Detecting d-wave superfluid and d-density wave states of ultracold Fermions on optical lattices 
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We propose a pump probe experiment for detecting the d-wave superfluid and d-density wave phases of ul- 
tracold Fermions on an optical lattice. The pump consists of periodic modulations of the optical lattice intensity 
which creates quasiparticle pairs in these systems. The changes in the momentum distribution under the drive 
can be used to measure quasiparticle dispersion and gap anisotropy. Further, we show that the pattern of peaks 
and dips in the spin selective density-density correlation function provides a phase sensitive probe of the sym- 
metry of the order parameter in these systems. 



Ultracold atoms on optical lattices are emerging as candi- 
dates for building a quantum simulator (TJ of various lattice 
models used to study strongly interacting many body systems. 
These systems have already been used to implement important 
lattice Hamiltonians like the Bose Hubbard model |2], which 
is central to the theory of superfluid-insulator transitions 
and attractive and repulsive Fermi Hubbard models, which 
are central to the study of strongly interacting superfluids (4) 
and high temperature superconductors ifTTIl respectively. The 
recent observation of Mott insulating states in the repulsive 
Fermi Hubbard model [6) shows a novel manifestation of 
strong interaction in these systems. 

The attractive feature of cold atom systems is that they pro- 
vide a clean setting with controlled and tunable access to the 
strongly interacting parameter regime of model Hamiltonians. 
However, a good quantum simulator also requires methods to 
extract information about the many body correlations in the 
system. This is specially relevant for simulating strongly cor- 
related systems with complex phase diagrams, high T c super- 
conductors, where the quantum simulation can shed light on 
the nature of the various phases and test the various compet- 
ing theories that have been proposed. A quantum simulator 
which aims to choose between competing theoretical descrip- 
tions should be able to probe the characteristic correlations 
and excitations of the various observed and proposed phases. 

In the setting of strongly correlated system condensed mat- 
ter setting like high T c superconductors, measurement of ther- 
modynamic quantities like specific heat and of transport prop- 
erties like resistivity, Hall conductivity, and Nerst effect as 
well as various spectroscopic techniques like neutron scatter- 
ing and optical conductivity have been used to identify the var- 
ious thermodynamic phases and have yielded a lot of informa- 
tion about correlations and excitations. Detailed information 
about the electronic structure obtained from ARPES ifTH and 
STM (HI Q2) measurements has been crucial in identifying 
the nature of these thermodynamic phases. A crucial aspect 
of constructing the phase diagram is to identify the non-trivial 
symmetry of the various order parameters e.g. the d-wave 
symmetry of the superconductor. The techniques mentioned 
above contain indirect information about the symmetry, but 
a more reliable method for probing the symmetry of the su- 
perconducting order parameter has been to perform a phase- 
sensitive measurement, i.e. look for a tt phase-shift in the cor- 
ner junction SQUID [ 13 ]. However, this method is not directly 



applicable to the cold atom setting. Hence, it is important to 
develop new experimental techniques for cold atom systems 
to probe the correlations and order parameter symmetry. 

In this Letter, we consider the problem of characterizing 
d-wave Fermionic superfluids and states with d-density wave 
order. The d-density wave state has been proposed as a possi- 
ble candidate for the pseudogap phase (9). We focus on two 
issues: (1) determining the quasi-particle dispersion and iden- 
tifying the location of the nodes in the gap function if there 
are any and (2) experimentally identifying the order parame- 
ter symmetry. Our method consists of looking at the response 
of the system to periodic modulation of the optical lattice in- 
tensity. This drive creates pairs of quasi-particles of opposite 
momenta predominantly with total energy corresponding to 
the frequency of the drive. The change in momentum dis- 
tribution rik in response to the drive can be used to measure 
the quasiparticle dispersion. Further, the pattern of peaks and 
dips in different density-density correlation functions, which 
result from the interference of the quasiparticles excited by 
the drive, encode the symmetry of the order parameter (e.g. 
s-wave vs. d-wave pairing). The momentum distribution can 
be measured by a ballistic expansion experiment, while the 
correlators can be measured by noise correlations in the real 
space density profiles. The Letter consists of two parts: in 
the first part, we describe our method and apply it to the d- 
wave superfluid case, and in the second part we apply it to the 
d-density wave case. 

d-wave Superfluid — We will use the simplified BCS de- 
scription of a d-wave superfluid, which helps in demonstration 
of our proposed methods; however the qualitative results that 
we obtain are independent of the microscopic description, and 
the non-BCS nature of the superfluid will only result in quan- 
titative discrepancies. Thus this method may be applied to the 
superfluid phase of the Fermi-Hubbard model (assuming that 
it exists). We start with the tight binding BCS Hamiltonian 
on a two dimensional square lattice, suitably generalized for 
arbitrary pairing symmetry, 

where c\. a is the creation operator for an atom with flavor 
a G {T, |} and momentum k and £k = — p is the 
tight binding dispersion with p the chemical potential and 
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FIG. 1: Brillouin zone of a tight binding d-wave BCS superfluid. 
Solid black lines indicate the location of the Fermi surface (with- 
out the pairing interaction), while the dotted black lines indicate the 
location of the nodes. The shading indicates the value of order pa- 
rameter Ak- Several equal quasiparticle energy contours, bananas, 
are shown with colored lines. 



e k -2J(cos 

kx H~ cos ky ) where J is the tunneling matrix 
element. The mean field Ak describes the pairing interaction, 
and its k dependence reflect the pairing symmetry. We focus 
on the case of d-wave pairing Ak = Ao(cos k x — cos k y ). 

We drive the system by periodic modulation of the optical 
lattice intensity V(t) = Vb + SV sin(o;t), which modulates 
the tunneling matrix element J — > J(l + gi (£)). The modula- 
tion of the superconducting gap Ao — > Ao (1 +#2 (t)) depends 
on the microscopic mechanism of pairing. However, the pre- 
cise form of g2(t) is not important for qualitative results; the 
important point is that both g\ (t) and #2 (t) have similar time 
dependence (~ s'm(ujt)) ifTTl . (For the purposes of the exam- 
ple calculations we set #2^) = 2gi(t) = 2Asin(cjt), where 
A ~ y/Vo/E r SVo/V. This choice is motivated by the expec- 
tation that A - J 2 /U ELO.) 

For illustrative purposes, we shall consider driving a sys- 
tem that starts at zero temperature, i.e. with the initial con- 
dition that all quasiparticle states are empty. As the pertur- 
bation Hamiltonian only creates or destroys quasiparticles in 
pairs, we truncate the Hilbert space to exclude states with un- 
paired quasiparticles. This truncated Hilbert space is spanned 
by wavefunctions of the BCS form 



*w = n (S^ + ^w4t c -ki) io> . 



(2) 



The initial state corresponds to the BCS solution ^k = 
sign(A k ) [\ ( 1+Ck/ff k)] 1/2 and t; k = [\ (1 - £ k /£ k )] 1/2 
with Ek = y^t;^ + A^ the quasiparticle dispersion of the un- 
perturbed state. The equations of motion for the coherence 
factors are 



id t 



A k (t) 
A k (t) 2(£ k + < 7l (t)e k ) 



(3) 



where A^(t) = A k (l + 92(f))- In describing the evolution, 
we do not maintain the self-consistency equations for A and 
ji. This is justified in the linear response regime where we do 
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FIG. 2: Change in the momentum distribution function under the 
influence of the drive (rik(t — 25/ J)) — (nfc(0)) for a BCS d-wave 
superfluid, with Ao = 0.5 J, for various drive frequencies: uo = 
A /2 (left), A (center), 1.5A (right). 



not create a lot of quasiparticles. 

To gain a better understanding of the dynamics, we use sec- 
ond order perturbation theory in A to study the evolution of the 
system. It is more convenient to perform perturbation theory 
calculations in the quasiparticle basis 

*(*) = II (V-LkW + ^2,k(t)7 kT 7l ki ) |*bcs), (4) 



where 7^ a — u\^c\. G + av\^c_ k - are the quasiparticle cre- 
ation operators. In the tight binding d-wave superfluid, the 
contours of equal quasiparticle energy look like bananas (See 
Fig [I]). When k is off resonance (i.e. away from the ba- 
nana contour corresponding to uo = 2£?k), with detuning 
o)k = 2E\t—ou, |^2, k| oscillate between and ~ Xu^v^e^/cu^ 
at approximately the drive frequency. On the other hand, if k 
is on resonance, then |V>2,k| oscillate between and 1 with 
the period Pk ~ (2A^k^kCk) -1 . Thus, after the system has 
been driven for some time, the quasi-particle momentum dis- 
tribution function (7 k 7k) w ih develop peaks along the banana 
contours corresponding to uj = 2E^, with maxima near the 
banana tips. To optimize the peak height, the system should 
be driven for a time t ~ P^ 1 for k near a banana tip. For the 
remainder of the paper, we shall solve the equations of motion 
numerically, but use the insights from this simple perturbation 
theory to qualitatively explain the results. 

Momentum distribution function nk — To extract infor- 
mation about the dispersion relation, we propose to look at 
the momentum distribution of the atoms, rik = (c k c k ) = 
2|^k(£)| 2 - Experimentally, rik may be obtained from a band- 
mapping free expansion experiment lfT8l . As an illustra- 
tion, we consider the case of a d-wave superfluid (/i = 4.5, 
Ao = 1.0, J = 1) that we drive at three different fre- 
quencies uj = {0.5,1.0,1.5}A . The change in momentum 
distribution at a time t = 25 is plotted in Fig [2] As ex- 
plained before, at a given frequency uj, the excitations are be- 
ing driven predominantly on the banana- shaped contour cor- 
responding to uj = 2E\z. This leads to large changes in the 
momentum distribution along the same contour. The disper- 
sion can then be extracted by varying the drive frequency. If 
the frequency decreases (increases), the banana- shaped con- 
tour shrinks (grows) as can be seen in the left (right) panel of 
Fig. [2] The fringes correspond to the diffraction pattern that is 
formed at short times. The dispersion will reveal the k depen- 
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FIG. 3: Left: Quasiparticle density \^, 2 (t = 25/J)| 2 . Red (blue) 
arrows indicate wave-vectors that connect tips of bananas between 
which Ak does not change (changes) sign. Right: Correlation func- 
tion {p^(q)pi(—q))t — {p^(q)pi(—q))o for a d-wave superfluid. 
Pairing symmetry is encoded in peaks (red) and dips (blue) that de- 
velop in the correlation function as described in the text, (for both 
panels to = Ao) 



dent gap function and hence the d-wave nature of the super- 
fluid. In this respect, this method is similar to the ARPES de- 
termination of the momentum dependence of the gap in high 
T c superconductors lfT9l . 

In this context we will briefly comment on the response of 
a s-wave superfluid to the same drive. Since the s-wave sys- 
tem is gapped, one would not observe any response before the 
drive frequency exceeds twice the gap. In contrast, the low 
energy nodal quasiparticles in a d-wave superfluid lead to re- 
sponse at all drive frequencies. 

Spin up-spin down density correlation function — Next, 
we show that the changes in the spin up-spin down density 
correlation function N ^ (q) = (p^(q)pi(—q)) show a pattern 
of peaks and dips which reveals the order parameter symme- 
try in a way similar to corner junction SQUID measurements 
in condensed matter systems. This correlation function can 
either be measured using elastic light scattering or extracted 
from the noise in real space measurement of the particle den- 
sity. The point of measuring this particular correlation func- 
tion is that it is sensitive to the symmetry of the pairing inter- 
actions in the system. 

We can relate N^(q) to the BCS wavefunction at time t 
defined by the coherence u\^(t) and v\^(t) as follows 

N n (l) = E (4 +qT W c kT(04,_ qi (i)c kU (t)} (5) 

k,k' 

= E (4 +qT W ct -k-<uW) (c-k|(i)c kT (i)> (6) 

k 

= £M*X(*K +q (*)W*)- (7) 

k 

Here, in the middle step, we have used the fact that for the 
BCS wavefunction only the anomalous pair correlations are 
non-zero. We note that the d-wave symmetry of u^v^. is pre- 
served throughout the time evolution of the system. The form 
of N^(q) looks like the product of Cooper pair wavefunc- 
tions at k and k + q, and therefore measures the interference 



of quasiparticles created at the respective wavevectors. Driv- 
ing creates predominantly creates quasiparticles near the tips 
of the resonant banana contours. The interference of these 
quasiparticles result in peaks or dips in the noise correlation 
function whenever the wave- vector q that appears in N^(q) 
connects a pair of banana tips located at k and k + q. 

The change of N^(q) after driving the system for ~ 4 cy- 
cles is depicted in the right panel of Fig. [3] for the case of the 
d-wave superfluid. To decipher the meaning of the various 
peaks and dips we compare the right panel of Fig.[3]to the left 
panel which depicts the corresponding change in the quasipar- 
ticle distribution function (7^ 7k)- Connecting the tips of the 
bananas (or the peaks in the quasiparticle density) with arrows 
in all possible ways we identify 7 unique wave-vectors, which 
are labeled q\ through qj in Fig. [3] The vectors that connect 
the remaining pairs of tips may be obtained by reflections of 
q\..qj across the nodal directions. The presence of a peak/dip 
at a given q vector depends on the time of drive; however the 
relative pattern of peaks and dips is governed by the symme- 
try of the order parameter. If the wave vector connecting two 
banana tips where Ak and Ak+ q has the same sign (say qi) 
has a peak in the signal, the wave- vectors connecting the tips 
where the order parameter has opposite signs (say c\2) should 
have a dip and vice versa. Note that for s-wave pairing where 
order parameter does not change sign we would have all peaks 
or all dips depending on time of drive. For other singlet pair- 
ing symmetries the relative pattern would change depending 
on the location in the Brillouin zone where the order parame- 
ter changes sign. 

Therefore, given the noise correlation function we can (1) 
determine the dispersion relation along the Fermi surface by 
measuring how the distances between the peaks and the dips 
varies with the drive frequency; (2) determine the pairing sym- 
metry by figuring out the relative sign of Ak at the various 
banana tips from the pattern of dips and peaks. 

The important point to note is that our method asks a bi- 
nary question of whether there is a peak or a dip much like 
corner junction SQUID experiment. In this sense this method 
is robust to quantitative changes in the signal due to strong 
correlation effects. 

d-Density Wave — Next, we show that similar methods can 
be applied to the d-density wave states which have been pro- 
posed as prospective candidates for the pseudogap phase of 
high T c superconductors [9 ). In this state the operator cj c c k+q 
acquires an expectation value Ak. We concentrate on the 
commensurate case in which the ordering occurs at momen- 
tum q = (7r,7r). At the mean field level the system is de- 
scribed by the Hamiltonian 

H o = Yl ^ kC k c k + A kc k c k+q + A k c k c k+q , (8) 

k 

where the mean field has the same form as for the d-wave 
superfluid case Ak = Ao(cos(k x ) — cos(k y )). The com- 
mensurate ordering halves the Brillouin zone but fills it with 
two different (quasiparticle) bands. The two bands touch at 
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FIG. 4: Effects of drive at frequency to = Ao/2 on den- 
sity wave states. Change in the momentum distribution func- 
tion (nk(t = 25/ J)) — (rik(O)) for d-density-wave at half- 
filling (top). Change in the density-density correlation function 
(/>kp q -k(t = 25/ J)) — (pkPq-k(t = 0)) for anisotropic s-density- 
wave Ak = 0.5| cos(fcc) — cos(k y )\ (left) and d-density-wave 
Ak = 0.5cos(fec) — cos(k y ) (right). 



the points (±tt/2, ±tt/2). At half-filling (/i = 0), the lower 
band is completely filled and the upper band is completely 
empty. Away from half filing there are empty states in the 
lower band. Optical lattice modulation excites quasiparticles 
from the lower branch to the upper branch leaving quasiholes 
behind in the lower band. In other words, in analogy with the 
d-wave superfluid, the drive creates quasihole-quasiparticle 
pairs with zero total momentum and total energy uo. The wave 
function in analogy to Eq. ([2]), can be written in the form 

*W = n( a kW4+/?kWc3 {+q ). (9) 

k 

The dispersion relation for quasiparticles can again be ob- 
tained by looking at the response of the momentum distribu- 
tion to the drive. Because of the folding, the momentum dis- 
tribution for the inner half-zone and outer half zone must be 
patched together 



n k 



M*)l 2 
IAc(t)| 2 



k G inner zone 
k G outer zone 



(10) 



There is one important difference in that tuning the system 
away from half-filling does not shift the location of the ba- 
nanas as the density wavevector q has been assumed to be 
fixed. So independent of filling, the centers of the bananas re- 
main fixed at (±7r/2, ±7t/2). Away from half-filling the low 
frequency response (inner portion of the banana) disappears as 
both branches are either empty or filled and thus a quasihole- 
quasiparticle pair cannot be created. 

The symmetry of the order parameter is encoded in the pat- 



tern of peaks and dips of the (pkPq-k)- We compare this 
pattern for two cases Ak = Ao(cos(/^) — cos(fc 2/ )) and 
Ak = Ao| cos(fcaj)— cos(fc 2/ )|. Because the peaks in the quasi- 
particle density near the tips of the bananas are not as sharp 
as for the case of d-wave superfluid it is hard to identify the 
origin of various peaks and dips, but there is a clear distinction 
between the two cases. 

Concluding remarks — Probing the properties of the corre- 
lated state of ultracold atoms is an important part of building a 
"quantum simulator." In this Letter we propose to use periodic 
modulation of the optical lattice amplitude to probe Fermionic 
superfluids and density wave states. This probe is relevant to 
simulating the repulsive Hubbard model which is thought to 
be a minimal model of high T c superconductors. The methods 
we suggest are useful for probing states with nodes or at least 
an anisotropic gap. We have shown that: First, we can deter- 
mine the quasi-particle dispersion and gap anisotropy by look- 
ing at the change in the momentum distribution created by the 
drive. Second, by measuring the appropriate density-density 
correlation function we can determine the pairing symmetry 
from the pattern of peaks and dips that develops. By mea- 
suring the distances between the peaks and dips, we can also 
extract the gap anisotropy along the Fermi- surface from the 
second method method. 
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